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В работе представлена нелинейная динамическая модель краткосрочного прогнозирования цен на фьючерсном рынке, 
разработанная на основе методов нелинейной динамики. Приведены схемы построения точечного и интервального прогнозов. 


Экономическая постановка задачи 

В данной работе объектом исследования явля- 
ется фондовый рынок. Динамика цены на фьюче- 
рсном рынке часто носит флуктуационный харак- 
тер, поэтому для описания данных процессов ис- 
пользуют стохастические вероятностные модели, в 
которых исследуемый процесс есть решение систе- 
мы стохастических уравнений, содержащих источ- 
ник случайности. 

Наиболее перспективным является метод, ос- 
нованный на теории детерминированного хаоса. 
Здесь возникновение флуктуаций объясняется как 
результат неслучайных взаимодействий связанных 
переменных в нелинейной динамической системе. 
Согласно данной теории введение в модель теоре- 
тически оправданных нелинейностей может опи- 
сать экономические флуктуации более успешно, 
чем введение случайных переменных [1]. 

На первом этапе моделирования необходимо 
определить основные факторы движения рынка. 
Согласно теории технического анализа, которая 
широко применяется для прогнозирования пове- 
дения рыночных характеристик, динамика рынка 
включает три основных источника информации, а 
именно: цену контракта, объем торгов и "открытый 
интерес". Объем торгов и "открытый интерес" не 
первостепенные, но, тем не менее, чрезвычайно 
важные факторы, влияющие на формирование це- 
ны. "Открытый интерес" - это количество не зак- 
рытых позиций на конец торгового дня [2]. На ос- 
нове данных факторов формируется ликвидность 
рынка и оборот торгов. 

В результате модель прогнозирования цен на 
фьючерсном рынке должна описывать процесс из- 
менения трех рыночных характеристик - цены 
контракта, объема торгов и "открытого интереса". 

Математическое обоснование и построение модели 

Экспериментально поведение сложной систе- 
мы определяется путем наблюдений в течение ка- 
кого-то интервала времени над некоторым эконо- 
мическим показателем Х(1), в нашем случае - це- 
ной. Анализ этой последовательности, на форми- 
рование которой оказывают влияние и другие пе- 
ременные, позволяет определить число дифферен- 
циальных уравнений первого порядка А г , необхо- 


димых для моделирования динамики системы. 
Фрактальная размерность аттрактора ё должна 
удовлетворять неравенству ё<К Округлив ё до 
ближайшего целого сверху, получим величину N 
[ 3 ]. 

Для определения размерности аттрактора стро- 
им псевдофазовое пространство, используя значе- 
ния временного ряда цены, взятые со сдвигом во 
времени. Например, фазовый портрет на плоскос- 
ти может быть построен с использованием вектора: 

ДО),ДН7)}. 

Идея заключается в том, что сигнал Х(і+Т) свя- 
зан с производной сигнала Х(1), и результат имеет 
те же свойства, что и при использовании истиной 
фазовой плоскости. 

Далее для численной оценки корреляционной 
размерности используем корреляционную функ- 
цию, подсчитывающую число пар точек, расстоя- 
ние между которыми меньше Ь [3]. 

С(Ы= Вт{Аг2еЦ--К-*ф>> 

ІѴ— »со ДГ И 7 11 

1 у іф у 

<2 = 0, если Ь - ||х - X / 1 > 0; 

0=1, если Ь - |х. - X ] I < 0. 

Сама корреляционная размерность оценивает- 
ся из наклона зависимости 1п[С(С)] от 1п(Х). 

Проведенный анализ последовательности наб- 
людений изменения цены на различных фондовых 
рынках позволяет ограничиться тремя уравнения- 
ми для описания исследуемой динамической сис- 
темы [3]. 

В качестве первой фазовой координаты выби- 
раем Хфі) - цену контракта, в качестве второй и 
третьей - рыночные характеристики, которые ока- 
зывают наиболее сильное влияние на формирова- 
ние цены - это Х 2 (1) - объем торгов и Хфі) - "откры- 
тый интерес". 

Поскольку исследуемые динамические процес- 
сы, в общем случае, описываются дифференциаль- 
ными уравнениями турбулентного типа, то модель 
должна включать систему трех нелинейных диффе- 
ренциальных уравнений следующего вида: 
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г 

= а, (I) ■ Х,(І) + а 2 (1) ■ Х,(І) ■ Х 2 (1) + а 3 (1) ■ Х,(1) ■ Х 3 ( 1 ) 
аі 

{ - ^ = Ъ,(і) ■ Х 2 (і) ■ Х,(і) + Ь 2 (І) ■ Х 2 (і) + Ъз(і) ■ Х 2 (і) ■ Хз(1) ( 1 ) 

аі 

ѢМ1= Сі(і ) . х,(1) Х,(1) + с 2 (1) ■ Хз(І) ■ Х 2 ( 1 ) + с 3 «) ■ Х 3 ( 1 ), 

^ аі 

Существующая взаимосвязь между вышеопи- 
санными экономическими показателями отражена 
перекрестным произведением соответствующих 
фазовых переменных: Х 1 (і)-Х 2 (1) - оборот торгов. 
Отражает взаимосвязь между ценой контракта и 
объемом торгов и позволяет учитывать в модели 
внутренние силы, управляющие движением цены. 

Х,((УХ 3 (() - текущая ликвидность рынка. Позво- 
ляет отразить факт заинтересованности тем или 
иным контрактом с долгосрочной точки зрения; 
другими словами, определить, насколько серьезно 
участники рынка воспринимают текущий тренд. 
Отражает взаимосвязь между ценой контракта и 
"открытым интересом". 

Х 2 (1)-Х 3 (Г) - взаимосвязь между объемом торгов и 
"открытым интересом". О количественной харак- 
теристике взаимосвязи объема и "открытого инте- 
реса" известно мало, но качественную характерис- 
тику, опираясь на экспериментальные данные, 
можно сформулировать следующим образом: "Уве- 
личение объема торгов должно подтверждаться 
достаточным открытым интересом" [2]. 

й/0, а 2 ((), я/ О, й/0, й/0, й/0, с/0, с/ О, с/ О - 
неизвестные коэффициенты системы, определяю- 
щие степень влияния соответствующих показате- 
лей рынка и их взаимосвязи на поведение системы. 
Данные коэффициенты являются переменными на 
некотором достаточно большом отрезке времени, 
но кусочно-постоянными на небольшом исследуе- 
мом интервале - шаге прогноза. Эти коэффициен- 
ты характеризуют: 

я/0 - изменение цены на один процент в еди- 
ницу времени, 1/с. 

й/0 - изменение оборота торгов на один про- 
цент в единицу времени, 1/с. 

й 3 (0 - изменение ликвидности рынка на один 
процент в единицу времени, 1/с. 

й/0 - влияние цены на изменение оборота тор- 
гов, 1/руб-с. 

й/0 - изменение объема торгов на один про- 
цент в единицу времени, 1/с. 

й/0 - влияние открытого интереса на измене- 
ние взаимосвязи: объем - "открытый интерес", 
1/шт-с. 

с/0 - влияние цены на изменение ликвидности 
рынка, 1/руб-с. 

с/0 - отражает влияние объема торгов на изме- 
нение взаимосвязи: объем - "открытый интерес", 
1/шт-с. 

с/0 - изменение открытого интереса на один 
процент в единицу времени, 1/с. 


Вычисление коэффициентов производиться по 
всем параметрам модели в фиксированные момен- 
ты времени. В результате получается система алгеб- 
раических уравнений (2) относительно неопреде- 
ленных коэффициентов. Первые производные в 
левых частях уравнений оцениваются при помощи 
кубического сплайна. 


^ ѢАА = а , . ед + а, ■ ед - ед + а, ■ ед ■ ед 
А, 


СІІ ( 

<ЬАЬ) _ 
д.1 { 

/ЧЙы) = 

Я м 

Я 2 0 м ) . 
А,., 


=ь,- ед ■ ед + ъ 2 ■ ед + ь 3 ■ ед ■ ед 
а -ед ед + с 2 ■ ед • ед + с, ■ ед 

- агХ/(Ь+]) + а 2 Х 3 ((і+і) Х 2 (Ь+і) + ауХ 2 (Ь+і ) -Х 3 (Ь+і) 

= Ь,Х 2 (и +І ) Х,« м ) + Ь 2 -Х 2 (1 М ) + Ъз-Х 2 (і м ) Х 3 (І М ) (2) 
= сгХ 3 (Ь+і ) Хі(и+і) + с 2 -Х 3 (Іі+і) ■ Х 2 (іі+і) + суХЛи і і) 


1 1+2 ^ — ауХі(Ь+ 2 ) ~ а 2 ш Х/(Ь+ 2 ) Х 2 (1 ІІ2 ) + а 3 Х;(1у 2 ) Х 3 (Ь+ 2 ) 

А,., 


= ь,-х 2 (Ьу 2 ) -х,(і 1+2 ) + ъ 2 х 2 (и +2 ) + ЪуХ 2 (и +2 ) -ед+д 

\ ~ СуХ 3 (іі+ 2 ) Х](іі+ 2 ) + С 2 Х 3 (Ь+ 2 ) Х 2 (Іі +2 ) + С;уХі(у 2 ) 

V А м 


А і+ 


Из решения данной системы находим искомые 
параметры, которые считаем постоянными на шаге 
прогноза. Подставляем а ь а 2 , а 3 , Ь и Ь 2 , Ь 3 , с,, с 2 , с 3 в 
систему уравнений (1) и, решая задачу Коши для 
системы обыкновенных дифференциальных урав- 
нений при начальных условиях в точке (?,+/, нахо- 
дим вектор прогностических значений Х р (х 3 , х 2 , х 3 ). 
В результате получаем точечный прогноз на один 
шаг вперед. 

Следует отметить, что все процессы, характери- 
зующиеся наличием хаоса, гиперчувствительны к 
точности задания параметров и начальных условий 
[1]. Поэтому краткосрочное прогнозирование наи- 
более качественно осуществляется с помощью 
адаптивных непрерывно подстраиваемых моделей. 
Это означает, что при прогнозировании с помощью 
модели (1) на каждом шаге осуществляется обнов- 
ление коэффициентов а,, а 2 , а 3 , Ь и Ъ 2 , Ъ 3 , с ь с 2 , с 3 и 
начальных условий с учетом развития событий. 


Построение интервального прогноза 

Предсказание вектора Х р одним значением на- 
зывается точечным прогнозом. Для построения 
интервального прогноза (предсказания того, что 
событие осуществиться в заданном интервале зна- 
чений с указанной вероятностью) необходимо, в 
первую очередь, проанализировать относительную 
ошибку прогноза, которая представляется после- 
довательностью - (е„ ( = 1, ... п). е, - должны быть 
случайными величинами и удовлетворять нор- 
мальному закону распределения. 

Анализируя временной ряд е„ доказываем слу- 
чайность колебаний уровней е„ пользуясь критери- 
ем пиков. Критерием случайности ряда в, с 5 % 
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уровнем значимости является выполнение следую- 
щего неравенства [4]: 


Р> [~(п~ 2) -1,96 


16/7-29 


90 


( 3 ) 


где Р - число поворотных точек ряда в,. Квадрат- 
ные скобки означают, что берется только целая 
часть полученного значения. 


т = Х Р ~ х йот (и + от - 2) 

(п - Р)Р> Ч , + (от - 1 )75 ѵ п + от 

где х р - предсказанные значения; х - реальные дан- 
ные; Б хр - дисперсия предсказанных значений це- 
ны, А - дисперсия реальных данных; и, от - соот- 
ветствующие относительные частоты, заданные 
статистической таблицей [5]. 


При выполнении данного критерия гипотеза о 
случайности колебаний е, принимается. 

Далее предполагаем, что е, имеет нормальное 
распределение, и рассматриваем вопрос о согласо- 
ванности статистического и теоретического расп- 
ределений, используя критерий Пирсона [5]. 


Величина Р характеризует вероятность согласо- 
ванности теоретического и статистического расп- 
ределений. При Р> 0, 1 нормальный закон достаточ- 
но удовлетворительно воспроизводит заданное ста- 
тистическое распределение. Р определяется из 
справочных таблиц по значениям % 2 и г [5]. 

г=1- 3 , 


Г 


к у(°Р Р,У 

Ы Рі 


(4) 


где щ - относительные частоты, заданные статис- 
тической таблицей; р, - вероятности, полученные 
по некоторому теоретическому закону распределе- 
ния, в нашем случае нормальному; / -число разря- 
дов статистической таблицы. 

В результате интервальный прогноз строится по 
следующей схеме: 


Ѵ Р = Х Р ± 



число и„ находят по справочным таблицам значе- 
ний функции Лапласа [5] из условия 


Ф(и с ) 


1 -а 

~1Г’ 


где ( \-а ) - доверительная вероятность; Х р - предс- 
казанное значение; А - дисперсия прогноза, кото- 
рая определяется по временному ряду г,. 


Проверка на адекватность и тестирование модели 

В общем случае под адекватностью понимают 
степень соответствия модели свойствам реального 
объекта или явления, для описания которых она 
строиться. 

В данной работе оценка адекватности модели 
проводилась по средним значениям откликов мо- 
дели и системы, то есть предсказанных значений 
показателей и реальных. При этом проверялась ги- 
потеза о близости среднего значения наблюдаемых 
переменных среднему значению реальных данных. 

По результатам испытаний вычисляем значение 
следующей величины [5]: 


Сравниваем Т„ с критическим значением Т„ 
взятым из справочной таблицы [5]. Если неравен- 
ство Т с > Т„ выполняется, то гипотеза об адекватнос- 
ти модели принимается с заданным уровнем зна- 
чимости. 


Таблица. Результаты прогноза цены контракта 


Истинное 

значение 
цены, $ 

Прогностическое 
значение цены, $ 

Относительная 
ошибка 
прогноза, % 

Доверительный 
интервал для 
прогностической 
цены 

171,50 

170,159 

0,781 

[170,137; 170,181] 

170,00 

169,460 

0,317 

[169,438;169,482] 

168,80 

168,328 

0,279 

[168,306:168,350] 

166,00 

165,458 

0,326 

[165,436; 165,480] 

163,70 

163,951 

0,153 

[163,929; 163,973] 

166,50 

164,671 

1,098 

[164,649; 164,693] 

158,40 

163,351 

3,215 

[163,329; 163,373] 

158,25 

160,714 

1,557 

[160,692;160,736] 

155,35 

160,514 

3,324 

[160,492:160,536] 

155,10 

157,156 

1,325 

[157,134; 157,178] 

158,20 

156,237 

1,240 

[156,215;156,259] 

155,05 

155,542 

0,317 

[155,520;155,564] 


Для проведения тестовых расчетов использова- 
лись исходные данные - итоги торгов с Чикагской 
товарно-сырьевой биржи по продажам кофе за нес- 
колько лет [6]. Результаты проведенного исследова- 
ния сведены в таблице и представлены на рисунке. 



Рисунок. Сравнительная диаграмма реальных и прогнозиру- 
емых 

Временной ряд в, для нашего прогноза удовлет- 
воряет приведенным выше критериям (3), (4) и 
имеет дисперсию Л р = 1 ,2829, поэтому нетрудно оп- 
ределить с вероятностью А=0,95, что относительная 
ошибка прогноза не превысит двух процентов. 

Р{е< 0,02) = 0,95 
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Для проверки адекватности модели по формуле 
(5) определяем величину Т„. Отмечаем, что 
неравенство Т с > Т„ выполняется, а это указывает на 
то, что модель адекватна описываемой системе с 
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